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Abstract 

The transition strengths for the four infrared-active vibrations of charged 
Cqo molecules are evaluated in self-consistent density functional theory using 
the local density approximation. The oscillator strengths for the second and 
fourth modes are strongly enhanced relative to the neutral Cqo molecule, in 
good agreement with the experimental observation of "giant resonances" for 
those two modes. Previous theory, based on a "charged phonon" model, 
predicted a quadratic dependence of the oscillator strength on doping, but 
this is not borne out in our calculations. 
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I. INTRODUCTION 



Recent experimental measurements |T]-^ of infrared absorption in A^Cqq (A = K, Rb), 
(x = 0, 1, 3, 4, 6) have demonstrated the strong dependence on doping of the oscillator 
strengths for several of the Ti„ infrared-active modes of the Ceo molecule. The second and 
fourth modes are enhanced by roughly 100, while the first and third show little or no change 
in oscillator strength in going from undoped Ceo to Ceo • 

The "charged phonon" model introduced by Rice and Choi obtains this increase 
in oscillator strength from the stronger coupling of the vibrational mode to an electronic 
transition. In particular the transition of the dopant electrons between the tiu and the 
higher-lying tig molecular orbitals could enhance the dipole associated with the vibrational 
modes. The result is, in addition to a softening of the phonon modes proportional to the 
number of dopant electrons x, an enhancement of the oscillator strengths proportional to x^. 
Friedman examined the effect of doping in a tight-binding model P,^. One parameter set 
did obtain the observed enhancement in the fourth mode, but the second mode could not 
be explained. Furthermore, experimental evidence for the x"^ doping dependence appears 
inconclusive. 

In a previous paper ^ we examined the oscillator strengths for the transitions in undoped 
Ceo, and found that the calculations using a frozen-phonon approach in the local density 
approximation gave good agreement with experiment on the magnitude of the oscillator 
strengths, particularly when compared with a tight-binding approach. In the present paper 
we apply the same self-consistent density functional technique to charged Ceo molecules, 
evaluating the dipole moment induced on the system by the same vibrational distortion 
used for the undoped system. 
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II. THEORY 



As described, we make an adiabatic approximation to calculate the oscillator strength 
associated with the vibrations. We will calculate by distorting the Ceo structure according to 
the atomic displacements in the vibrational modes, and then evaluating the dipole moment 
of the structure. Defining the displacement of an atom i in the mode a by the vector df, 
the oscillator strength is given by the formula 

McE^df■df 

where Da is the dipole moment in the distorted structure and mg, Mc are the electron and 
carbon masses, respectively. The displacement vectors are the same as we used for undoped 
Ceo 0' calculated from the phenomenological model of Ceo vibrations in ref. Equation |l| 
is of course the standard way to obtain absorption strengths, although it is usually presented 
in terms of infinitesimal displacements {i.e. the dipole derivatives), rather than the finite 
displacements we use. 

The electrons of the system are described by the Kohn-Sham density functional approach 
0, in which an effective one-electron problem is introduced that depends on a self-consistent 
charge density through an approximate "exchange- correlation" functional, which takes into 
account the electron-electron interactions. We use the local-density approximation (LDA) 
10| for this exchange-correlation functional, which has proven remarkably successful in de- 



scribing the charge distributions and energies of a wide variety of solid and molecular ma- 
terials. Our calculation of the electron charge density for these fixed ionic positions uses a 
version of the Car-Parrinello electronic structure technique implemented in a code 

provided by Wiggs |jl3|, recently modified to facilitate rapid convergence to the ground state 
The same pseudopotential [|15| and plane-wave basis set (with 35 Ry cutoff) was used 



here as in our earlier work. The discretization of the plane wave basis set is equivalent to 
treating the problem in a cubic box with periodic boundary conditions. The equivalent box 
length in our calculations was 23 atomic units, or 12.2 A. Wave functions for the charged 
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Ceo molecules were calculated self-consistently for the excess charge and distorted molecule 
conditions, with double occupancy assumed for all the lowest-energy orbitals, except that 
the highest occupied level contained only one electron for the odd-charge molecules. Note 
that the periodic boundary conditions imply that the Hamiltonian effectively has a uniform 
positive background charge, instead of the discrete positive ion charges that would be present 
in the real M^^Ceo material. 

This technique was then applied to calculate the charge density and dipole moment for 
each of the frozen Tiu modes for the Ceo ions, using a displacement of net amplitude of 
lA, which in our earlier work was shown to be still within the regime where the moment 
varies linearly with the displacements. We estimate the random errors (due to insufficient 
convergence or box size) in the resulting dipole moments to be roughly 0.1 eA, or between 
6% for the largest moments and 100% for the smallest {ie. the smallest /-strengths we 
calculate are not significantly different from zero). 

Before proceeding to the oscillator strengths, we note that the vibrational frequencies are 
readily obtained from the energy differences due to these displacements, and we have tabu- 
lated the results for neutral and charge -6 Ceo in Table |. The pronounced mode softening 
apparent in the experiments for the first and fourth modes is absent in our calculated results 
- in fact we find the first mode increases in frequency. This is most likely a consequence of 
our assumption of identical vibrational displacements in the neutral and charged systems. 

The results for the oscillator strengths are much more conclusive. As can be seen from 
Figure |l], none of the oscillator strengths we obtain increases quadratically with x. In 
agreement with experiment, the second and fourth modes are strongly enhanced in strength. 
Our calculations find the first and third to actually diminish in oscillator strength with 
doping. The experiments do not clearly demonstrate this, although the fact that these 
peaks are not observed in some of the experiments at intermediate doping is evidence that 
they may have decreased in amplitude. 

Figure ^ shows the relation between our calculations for the neutral and charged Ceo 
molecules and the available experiments on the undoped and doped systems. We have con- 

4 



verted the experimental measurements to units corresponding to our theoretical isolated Ceo 
/-strengths as outlined in the appendix. Theory and experiment are generally in agreement 
to within a factor of 2, which is remarkable because the strengths vary over 2 orders of mag- 
nitude, and different experiments themselves have discrepancies of roughly the same size. 
However, there are a few modes which were observed experimentally to have appreciable 
strength even though the theoretical strengths were much smaller - in particular the first 
mode for Cgo and the first and third modes for Cgo ■ This could be due to a breakdown in the 
theoretical approximations, or possibly due to contamination of the experimental samples 
by some fraction of neutral or charge -1 Ceo- The most significant disagreement appears to 
be for the fourth mode of the charge -6 species, where the theoretical strength is between 
3 and 5 times smaller than the experimental numbers. The most likely explanation of this 
is the breakdown of our theoretical assumption that the atomic motions in the vibrational 
modes are independent of charge state. However there may be other approximations (such 
as ignoring the actual positions of the compensating charges, or the conversion from solid 
to molecular strengths) that are invalid for this particular data point. 



III. CONCLUSION 

Our calculations based on the Kohn-Sham density functional theory thus appear to be 
in excellent agreement with the available experiments. In particular, this self- consistent 
approach correctly predicts the strong enhancement of the second and fourth modes in the 
doped systems. The quantitative agreement is generally within a factor of two, which is 
also roughly the extent to which different experiments agree. However, all but one of the 
modes for the Ceo significantly weaker in the theoretical calculations than in the 

experiments, possibly indicating that for such a highly charged ion the approximations we 
have made are no longer suitable. The experimental observation of the weak first and third 
modes could also be ascribed to sample inhomogeneities and resulting contamination from 
other charge states, however. 
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We note here that the net dipole moment for a particular oscillation mode results from 
a very delicate cancellation of contributions from almost all of the electronic states in the 
system - most of the states have contributions to the total dipole that are roughly of the 
same magnitude (one electron-A) as the final dipole moment. Thus the charged-phonon 
model of Rice and Choi [^] appears to be oversimplified. As shown in our Figure the 
quadratic doping dependence of the strengths predicted by that model is not observed for 
any of the modes in the theoretical calculations. Unfortunately, our calculations do not 
suggest any simple explanation of the doping dependence - it appears that the only way to 
obtain meaningful theoretical numbers is through such a self-consistent calculation of the 
response of the electrons to the vibrations. 
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APPENDIX: THE EXPERIMENTAL OSCILLATOR STRENGTHS 

In the published experimental work, the absolute strengths of the absorption peaks are 
generally expressed in terms of an effective plasma frequency uip for each mode, ie. a 
contribution 

A-KxM = ujI/{uj] -uj^- tUJ^j) (Al) 

associated with mode j to the dielectric function of the solid material, where Uj is the 
observed frequency of the resonance. Theoretically, the dielectric response e{uj) of the solid 
derives from the polarizability 7moi('-^) of ^^e single molecule through the Clausis-Mossotti 
relation: 
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e{u) = 1 + 47r- 5-— — — - (A2) 

1 - —J\^niol[i^) 

where N is the density of molecules in the solid. The fundamental /-strengths we cal- 
culate are associated with the single-molecule polarizability (ignoring damping) due to all 
electromagnetic resonance modes in the molecule: 

The denominator in equation (|A2|) shifts the resonant frequency from UjQ to cjj, and also 
changes the overall strength by a factor depending on the other modes. These other modes 
can be absorbed into a real net background contribution to the dielectric function that 
manifests itself simply as a value of the refractive index n different from 1 (recall e = 
{n + ikY). The resulting relation between the effective plasma frequency ujp for a mode and 
the associated molecular /-strength fj is: 

£ 9 TfX 2 / A /I \ 

^ (n2 + 2)2 47riVe2'^^- ^ > 

We have used equation ( [A4| ) directly to convert the results of reference 0. The ex- 
perimental results of reference requirement multiplication of the tabulated S by the 
vibrational frequency uJq to obtain uj'i, and then fj. Reference |]16| was treated in our pre- 



vious paper 0. Reference |jT| only provided relative strengths between the charge -6 and 
neutral species: we assumed the numbers for the neutral species were identical to that of 
reference p|, and plotted the resulting /-strengths for the charged molecule. 
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TABLES 



TABLE I. Energies of the vibrations, in eV. 



Label 


Ceo 




*-^eo 




Experiment 


LDA 


Experiment 


LDA 


Tiuil) 


0.065 


0.065 


0.058 


0.069 


Ti„(2) 


0.071 


0.070 


0.070 


0.070 




0.147 


0.149 


0.147 


0.146 


ri„(4) 


0.177 


0.197 


0.166 


0.195 
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FIGURES 

FIG. 1. /-strengths per Ceo from our LDA calculations for the neutral and charged molecules 
with up to 6 extra electrons, for the four Tiu vibrational modes. The curves are labeled by the 
associated vibrational mode. 

FIG. 2. Comparison of our calculated strengths with published experimental measurements 
for each of the vibrational modes for (a) neutral Cgo, and charges of (b) -1, (c) -3, (d) -4, and (e) 
-6. The experimental numbers come from references ||T^ (circles), (bursts), Q (diamonds), and 
IQ] (squares). 
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